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Abstract. We consider nonclassical entropy solutions to scalar conservation laws with concave- 
convex flux functions, whose set of left- and right-hand admissible states Ui,Ur across undercom- 
pressive shocks is selected by a kinetic function Ur = 'p^(ui). We introduce a new definition for the 
(generalized) strength of classical and nonclassical shocks, allowing us to propose a generalized notion 
of total variation functional. Relying only upon the natural assumption that the composite function 
(f^ o tp^ is uniformly contracting, we prove that the generalized total variation of front-tracking ap- 
proximations is non-increasing in time, and we conclude with the existence of nonclassical solutions 
to the initial-value problem. We also propose two definitions of generalized interaction potentials 
which are adapted to handle nonclassical entropy solutions and we investigate their monotonicity 
properties. In particular, we exhibit an interaction functional which is globally non-increasing along 
a splitting-merging interaction pattern. 
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1. Introduction. Consider the following initial- value problem associated with 
a conservation law in one-space variable : 

ut+f{u)^ = 0, ^^^^ 
u{0, •) = uo, 

where uq is a function of bounded variation on R, and the (smooth) flux / : M ^ R 
is a concave-convex function, in the sense that 

uf"{u)>0 {u^Q), /"'(0)^0, 

lim f'{u) = +00. (1-2) 

\u\ — *-t-oc 

Following LeFloch [B] , we consider nonclassical entropy solutions to this problem. 
Recall that, in many applications, only a single entropy inequality can be imposed on 
the solutions, i.e. 

U{u)t+F{u)^<Q, (f.3) 

where the so-called entropy t/ is a given, strictly convex function and the entropy flux 
F{u) :— J U'{v)f'{v)dv is determined by U. It is not difficult to construct multiple 
weak solutions to the same initial- value problem (|f .f p . (|f .3|) . so that one realizes that 
the single entropy inequality is too lax to determine a unique weak solution. In fact, 
for initial data restricted to lie in one region of concavity or convexity, the classical 
theory applies and leads to a unique entropy solution. Non-uniqueness arises when 
weak solutions contain transitions from positive to negative values, or vice- versa. 
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The above non-uniqueness property is closely related to the fact that discontin- 
uous solutions, in general, depend upon their regularization, that is, different regu- 
larizations or approximations to the conservation law (jl.ip may well lead to different 
solutions in the limit. This, in particular, is true for solutions to the Riemann problem, 
corresponding to the special initial data 



where ui,Ur are constant states. Indeed, for a wide class of regularizations, including 
regularizations by nonlinear diffusion-dispersion terms, nonclassical shocks violating 
Oleinik's entropy inequalities exist. The selection of a nonclassical solution is (essen- 
tially) equivalent to prescribing a kinetic function (/?'': M —s- R which, by definition, 
provides a characterization of admissible nonclassical shocks connecting states w_, w+, 
that is. 



For scalar conservation laws and, more generally, nonlinear hyperbolic systems, 
LeFloch and co-authors initiated the development of a theory of nonclassical entropy 
solutions selected by the kinetic relation (jl.Sp : various analytical and numerical as- 
pects have been covered. We refer to [6] for a review of the theory and to [8j for recent 
results on the numerical issues. On the other hand, the kinetic relation was originally 
introduced in the context of a hyperbolic-elliptic model describing the dynamics of 
phase transitions in liquids or solids, for which we refer the reader to Slemrod [12], 
Truskinovsky TS], and Abeyaratne and Knowles [T]. LeFloch [S] used the Glimm 
scheme and first investigated mathematically the kinetic relations arising in phase 
dynamics. The subject has developed extensively since then, and we will not try here 
to review the literature. 

Our objective in the present paper is to discuss the construction of new functionals 
for the (generalized) total variation and wave interaction potential of nonclassical en- 
tropy solutions to We consider solutions generated by Dafermos' front-tracking 
method, when the local Riemann solutions are nonclassical and are determined by a 
given kinetic relation (jl.Sp . We are interested in deriving uniform estimates for the 
total variation of solutions and showing that the scheme converges to global-in-time, 
nonclassical entropy solutions to the initial- value problem 

The new definition of generalized total variation proposed here is completely 
natural, and uses in a direct manner a key property satisfied by the kinetic function, 
that is, the uniform contraction property of the second iterate (p^ o ip^ . Indeed, this 
property (see (|2.ip . below) motivates our definitions of shock strength and, hence, 
of total variation. In turn, the proposed construction provides both an improvement 
of earlier results on the subject and a drastic simplification of earlier arguments (cf., 
for instance, with Section 8.1 in [6]). Most importantly, it appears that our new 
arguments are robust and may be generalized to tackle systems of equations. 

An outline of this paper follows. In Section [2l we begin with a brief review of the 
theory of kinetic relations, by emphasizing what we will need in the rest of this paper. 
We then introduce our new definition of generalized shock strength. In Section [3l we 
prove that the proposed generalized total variation functional is non-increasing along 
a sequence of Dafermos' front tracking solutions; cf. Theorem 13.11 In Section [H we 
turn to the construction of interaction functionals adapted to nonclassical entropy 




X <0, 
X > 0, 



(1.4) 



u+ = (p (u_). 



(1.5) 
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solutions. In Theorem 14.21 we show that a rather natural definition leads to an inter- 
action functional which is non-increasing in all but four interaction patterns. Next, in 
Section [5] we prove that the proposed interaction functional is actually globally non- 
increasing, at least in the significant case of merging-splitting wave patterns. Such 
solutions were originally introduced by LeFloch and Shearer [9]. Hence, Theorem 15. II 
below demonstrates the relevance of the proposed functional to handle nonclassical 
entropy solutions. Finally, in Section[6]we investigate a second definition of interaction 
functional which may be better suited for the extension to systems; see Theorem 14.21 

2. Kinetic function and generalized wave strength. 

2.1. Assumptions on the kinetic function. The general theory of kinetic 
functions if)' imposes the following conditions on the kinetic function: 

(Al) The map : M ^ K. is Lipschitz continuous and one-to-one. 

(A2) du'f^{u) < for all u, and ip^{0) = 0. 

(A3) The second iterate ip^ o is a strict contraction, i.e. for some K £ (0, 1) 

\ip^ oip\u)\ < K\u\, u^O. (2.1) 

These are the sole conditions required on the kinetic functions in the present paper. 
We recall that these assumption are satisfied by quite general regularizations of the 
conservation law (jl.ip and, moreover, natural analogues of the properties (Al), (A2), 
and (A3) are known to hold for systems, when ip^ is properly defined in this more 
general setting. We refer the reader to the monograph [5] and the references cited 
therein for further information. 

The kinetic function determines, via (|1.5[) . the set of admissible states 
arising on the left- and the right-hand sides of a discontinuity. When u_ > (m_ < 0, 
respectively), nonclassical shocks are required when the amplitude of the shock is 
large enough and the threshold u+ < ip^u^) {(p^u^) < u+, resp.) is reached. This 
threshold function lys" : R ^ R is defined as the unique value (p^{u^) ^ ip^{u^)^ 
such that 

f{u.)-f{^'{u^)) /(u_)-/(y«(^._)) 

U-—(p^{U-) (y9t'(u_) 

Geometrically, ((^"(m-), /((^"(m-))) is the point at which the line between /(u_)) 
and (^ip^{u-), f(^(p^{u-))) crosses the graph of / and therefore, if for instance U- > 0, 
the following inequalities must hold 

ip\u^) < cp^u^) <u^. (2.2) 

Under the concave-convex condition (|1.2p . the nonclassical entropy solution to 
the Riemann problem with data (|1.4p . > 0, is given by 

i) a shock if (p'^{ui) < Ur, or 

ii) a nonclassical shock connecting ui to (p^{ui) followed by a classical shock 
connecting (p^{ui) to Ur if ^^{ui) < Ur < (p^{ui), or 

iii) a nonclassical shock connecting ui to (p^{ui) followed by a classical rarefaction 
connecting 'p^{ui) to Ur if Ur < tp^{ui). 

When u; < and > 0, the nonclassical Riemann solver is similar; see [5] for further 
details. 
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Some remarks on our assumptions (Al)-(A3) above are in order. In fact, these 
properties can be motivated from similar properties satisfied by the so-called zero 
dissipation kinetic function ipQ, which is determined from the flux / and the given 
entropy U and is formally associated to a dispersion-only regularization of (jl.ip . 
Interestingly, this special kinetic function satisfies a stronger version of (j2.1|) , namely 

ip'o0^liu)=u. (2.3) 

On the other hand, for (non-trivial) regularizations including diffusion terms, the 
kinetic function satisfies the strict inequality 

|^^(u)| < |^^(u)|, u^O. (2.4) 

In the general theory, (j2.ip arises as a consequence of the more fundamental identities 
(|2.3p and (|2.4p . It is, therefore, natural that our functional involve the function i^jg 
rather than directly (p^ , and its definition makes direct use of the properties (|2.3p and 
(EH). 

Finally, we observe that for any bounded set of values \u\ < ill, properties (A2) 
and (A3) imply that there exist lower and upper Lipschitz constants Lip{u — ipQOip^) > 
and Lip(M — ip^^ o ip^) < 1 such that for u ^ 

Lip(ii — (fiQ o (p^) |m| < |u — (y9Q o (/?''(w)| < Lip(u — ifiQ o ip^)\u\. (2-5) 

In the context of the general theory, we note that only the lower bound is non- 
trivial and that it does hold for large classes of diffusion-dispersion relations. On the 
other hand, for highly degenerate regularizations, such as the one studied in [5], this 
condition might be violated. 

2.2. New notion of generalized wave strength. We introduce here a new 
functional, which is equivalent to the total variation and can be expressed entirely in 
terms of the kinetic function tp^ . In order to hope for an extension to systems, the 
formulation of the total variation should let the property (|2.ip appear as naturally 
as possible. Intuitively, the following definition of strength is based on the idea of 
mapping all the states to the convex region of the flux (u > 0), using ipQ. 

Definition 2.1 (Notion of generalized wave strength). To each classical or 
nonclassical wave (ii-,M+) one associates the generalized strength a{u-,u+) defined 
as follows : when u- > 0, 

\u--u+\, u+>0, ^2g) 

|u_ - (^^(w+)|, -U+ < 0, 



l'^o("-)-"+h > 0. 




The following notation will be helpful for the analysis. Henceforth, classical shocks 
and rarefactions joining two positive states is denoted by Cj_ and R\_, respectively. 
Similarly, when both neighboring states are negative, we write Cl and When a 
shock joins a positive state with a negative state, we write or depending on 
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whether or not the shock is classical or nonclassical, respectively. When the signs of 
the neighboring states are reversed, we simply write and . Shocks Cj_ and 
are also sometimes called crossing shocks. 

The generalized wave strength enjoys several immediate properties. 

• The first important property to observe is that the proposed generalized 
strength is continuous as u+ crosses ip'^{u-) and the solution of the Riemann 
problem goes from a single crossing shock to a nonclassical shock followed by 
a classical shock. 

For u- > 0, this is checked from the inequalities < ip^u^) < and 

the property (|2.1|) which yield 

In fact, (^"(u-) could be positive and a similar computation would show that 
continuity still holds. 

• The second important property of the generalized strength is its "equivalence" 
with the usual notion of strength. For positive non-crossing shocks and rar- 
efactions, the definition is the same as the usual one. When the rarefaction 
and the non-crossing shocks have two negative neighboring states, then there 
exists a positive constant C such that 

<j{u^,u+) = l^^oiu^) - > C>_ - u+\, (2.8) 

as long as stay within a bounded neighbourhood of the origin. For 

the shocks and iV^, it suffices to use (12.51) to show that the definition is 
equivalent to the usual notion of strength : 

a{Ni) = |u_ -(^^o(^^(m_)| > C"2\u- \ > C""|u_ -(^^(u_)|. (2.9) 

The same argument also applies to shocks and N^. Interestingly, our 
definition of shock strength makes clear the intuitive idea that it should be 
increasingly difficult to measure the strength of nonclassical shocks as the 
zero-diffusion (dispersion-only) case is approached. 

3. Diminishing total variation functional. We now introduce front-tracking 
approximate solutions to p.ip based on a nonclassical Riemann solver, as was con- 
structed by Dafermos ^ in the classical situation. These approximations are piecewise 
constant in space and are determined from the nonclassical Riemann solver described 
in the previous section. 

The first step of their construction is to build a piecewise constant approximation 
of the initial data uq which admits finitely many discontinuities and approaches uq in 
the norm with an error e, for some small e. The Rankine-Hugoniot condition can 
be used to propagate, in a conservative manner, the discontinuities of the initial data. 
When the Riemann solver calls for continuous waves, we replace them by a sequence 
of small discontinuities (u_,m+) whose strength satisfy a(u_,M+) < e. 

When two discontinuities meet, the nonclassical Riemann solver is used to deter- 
mine the neighboring states, but we continue to enforce that all outgoing waves be 
discontinuities. Despite the fact that the number of outgoing waves may be larger 
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than the number of incoming discontinuities, one can check (see [5] for details) that 
the total number of discontinuities remains bounded for all times, so that the front- 
tracking approximation can be defined for all times. 

For a front-tracking approximation u : M-|- x R M, formed entirely of propagat- 
ing discontinuities (each denoted by a), the inequalities (|2.8p and (|2.9|) . as well as the 
fact that the kinetic functions are Lipschitz continuous, imply that 

V{u{t)):^Y.^{u'^-^<)^ (3.1) 
is equivalent to the total variation norm 

TV(u(f)) := ^ (3.2) 

a 

where denote the left- and right-hand states of the discontinuity a. 

Incidently, we mention that the kinetic function (^q induces an isomctry on the 
space of BV(R) functions with respect to the norm V : 

a>^: BV(R) ^ BV(R) 

Using (|2.3p , it is straightforward to check that <i>Q is an isomctry. This result will not 
be used directly in this paper. 

Theorem 3.1 (Diminishing total variation functional). Let be a kinetic 
function satisfying the properties (A1)-(A3). For every front-tracking approxima- 
tion u : IR+ X M — > R to the conservation law (|l.ip based on the nonclassical Riemann 
solver associated with , the generalized total variation functional V(^u{t)'j is non- 
increasing. Precisely, the change in V during an interaction is given by: 

Cases RC-1, RC-3, CR-1, CR-2, CR-l 

Cases RC-2, RN, , , 

(3 4) 

cr(i?°"*)). Case CR-3, ^ ' ' 

all other cases. 

Here, ci = min{l, Lip((y5Q)}, C2 = Lip(M — i^q o ip^), while R™ and denote the 

incoming and outgoing rarefactions at an interaction, respectively. (The list of inter- 
actions is specified in the proof below.) 

From this theorem we immediately deduce the following: 

Corollary 3.2 (Existence of nonclassical entropy solutions). For any initial 
datauo G L°°(R)nBV(R) of (jl.ip . and any sequence of front-tracking approximations 
u'^ such that u^(-,0) converges to uq{-) in L^(R), there exists a subsequence of front- 
tracking approximations that converge in Lip ([0, T), L"'^(R)) n L°°([0, T), BV(M)) to a 
solution of the initial value problem (jl.ip . 

Proof of Theorem \3.1[ Our proof is a generalization of the one given in [6] (see 
Section 8.1). We need here to compute the variation of our functional V by distin- 
guishing between 16 possible interactions, after assuming ui > for definiteness, as 
described in Section 4.3 of [6 . We note that we have one more case to consider here, 
since our assumptions on the kinetic function are general enough to allow for the so- 
called CC-3 interactions. The subscript ' is used to indicate that a wave is outgoing 
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and with some abuse of notation for the wave strength we write, for example, a{N^) 
for the shock strength of a nonclassical wave. During a generic interaction between 
two waves, we denote the states on both sides of the left-hand wave by ui and Um 
while those associated with the right-hand wave are denoted by Um and Ur. Finally, 
the bounds on the change of V depend only on the properties (|2.ip - (|2.5[) which are 
therefore used freely throughout. 

Case RC-1 : {rIc^)-{C^'). This case is determined by the constraints 

mSL-K(ip\ui),(p\um)) <Ur < Ui, < < Um- 

This is further subdivided into two subcases depending on the sign of Ur- When 
Ur > 0, then the interactions are entirely classical {R\_C\_)-{C\^ ) and the 
inequalities < Ur < ui < Um suffice to check that 

[V] = - cj{r\) - a{Ci) = -2\u^ - ui\ = ~2a{Rl). 

When Mr < 0, then the interaction involves crossing shocks {R\ci)-{Ci ) and the 
states involved in measuring the strengths of the waves are 

< <p\){Ur) <Ui < Um, (3.5) 

since (p\^{ur) < tp\^ o ip^{ui) < ui. Then 

[V]=a{ci)-<r{Rl)-a{Ci) 

= \ui - (Po{Ur)\ - \u„i - Ul\ - \u,n - 'Poi'^r)\ = -2a-(i?^). 

Case RC-2 : {RlC^)-{Ni' R^'). This case is defined by 

ipHum) < Ur < (p^{ui) <0 < Ul < Um- 

In the first subcase, we assume (f\^{ur) < ui and use the previous conditions to 
deduce 

(PO O (p^{ui) < (poiUr) <Ui < Um- (3.6) 

The analysis when ui < Lp\^{ur) will not be treated since only the order of those two 
terms changes but the conclusions remain the same. 

[V]=a{Ni)+a{Ri)-a{Rl)-a{Ci) 

= \ui - ip\^0 ip^{ui)\ + |<^o(Ur) - V'O ° - |Wm - M/| - \Um ~ ^liUr)\ 

= 2|<y9o(Ur) - O ip^{ui)\ - 2\Um -ui\. 

Since (p'^{um) < Ur < 0, from property (|2.ip we deduce 

< ipQ{ur) < (fQ o tp^{um) < ^^0° f^{um)- Combining this with p.6p . we find 
[V] < -2Lip(u - (/Jq ° f^)\um -ui\ = -2Lip(u - (fQ o (^^)cr(i?+). 

Case RC-3 : {RlC^)-{Ni' C^'). The conditions initially satisfied by the 
neighboring states of the incoming waves are 

max {ip^{ui), ip^{Um)) <Ur < ip^{ui), < Ui < Um- 
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In the first subcase, we assume Ur <0 and the interaction is {R\_Cj.)-{N^ cL ) with 
the following states appearing in the strength of the waves 

iPoiUr) < ° <Ul < Um- (3.7) 

Therefore, the change in V is 

[V]^a{Ni)+a{Cl')-a{Rl)-a{Ci) 

= \UI -(PqO (p'iui)] + \(Po O ip^{ui) - ^Po{Ur)\ - \Um-Ul\ - |u„ - <^o("»-)l 

= -2\um-ui\ = -2a{Rl). 

In the second subcase with < Mr, we have Ur < f^iui) < ip^ o ip^{ui) < fySg o iy9^(u;) 
and therefore 

< <^o ° '^''("0 <Ul < Um- (3.8) 

The change in V is then 

[V] = a(Al') + a (4') - a{Rl) a{Cl) 

= \UI -tPoO (p\ui)\ + \(Po O (p\ui) - Ur\ -\Um-Ul\- \Um - "rl 

= -2\um - ui\ = -2a{Rl). 
Case RN : {RlNi)-{Ni R^ ). The states on both sides of the waves satisfy 

< Ul < Um and Ur = if^{Um)- 

Two cases again occur depending on the relative order of ui with respect to 
^g(Mr) = <^o ° V^(^m). We consider only the case where 

IPQ O ip^{ui) < (PqO ip^{Um) <Ui < Urn, (3.9) 

the other, (p^ o (p^{ui) < m; < tpg o ip^(um) < Um, being similar. In this case, the 
change in V is 

lV]=a{Ni)+a{R{_')-a{Rl)-a{Ni) 

= \UI - ipoO ip^{ui)\ + |(/?o ° 'P^iUm) - 'P^{.Ul)\ 

- \Um - Ui\ - \Um -'Po° 'p\um)\ 
= 2|<Po ° V'iUm) -^0° V'iui)] - 2\Um - Ul\ 

< — 2Lip(w — (PqO (p')\um — ui\ = — 2Lip(u — <^o ° V^)o'(-R+)- 

Case CR-1 : (Ciiii)-(Ci ). This is a simple case where the states are initially 
ordered as 

(p^{ui) < Ur < Um <0 < Ul. 

This provides ipQ{um) < foiur) < '■Poi'^H''^i)) < ''^i- The result is : 
[V]=aiCi')-<j{Ci)-a{R^) 

= Wl - foMl - Wl - (Po{Um)\ - IfaiUr) - <Po("m)l 
= -2\^\,{Ur) - ipl{Um)\ < -2Up{foMR-)- 
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Case CR-2 : {C\_r\_)-{C\_ ). The waves are entirely classical since 

< Um < Ur < ui. There is nothing new to check and the change in V is 

immediately found to be 

[V] = -2\ur-Um\ = -2a{Rl). 
Case CR-3 : {CiRi)-{Nl' Ri_'). The states begin in the order 

Ur < (f'iui) < ip^{ui) < Um < < Ul. 

One subcase is obtained when we assume that ui < (^o(^r)- The important states 
appearing in the definition of the strength of waves are then ordered as 

<^o(«m) < ° V^WO < < <^oK)- (3.10) 

The change in V is then found to be 

[V] = <j{Ni') + <j{r{_') - <j{ci) - 

= Wl - <^o ° V^Ml + \ foM ->fo° f^Ui)] - \ui - <fio{Um)\ - \<foiUr) " <fo{Um)\ 
= -2\ipl o (p^{ui) - (f>l{Um)\- 

In the second subcase, when <Po('Wr) < we use 

^o{Um) < </2o O ip''{ui) < ifoiUr) < Ul (3.11) 

to deduce the same equality : 

[V] = \UI -foO <f\ui)\ + \ fo{Ur) -IfQO ip\ui)\ - \UI - <fo{Um)\ " \<fo{Ur) - <fo{Um)\ 

We now observe that in both subcases, the change is equal to a physically relevant 
quantitiy 

[y] = -2(a(i?i)-a(i?i')). 
Case CR-4 : (C^i?i)-(7V^ Cl ). The states of the incoming waves satisfy 

(p^{ui) <Ur < (fi^Ui) < Um <0 < Ul. (3.12) 

This leads us to the inequalities 

^o{Um) < Vo("r) < </?0 ° ^H^l) < Ul- (3.13) 

The change in our functional is therefore 

[V]=a{Ni)+a{Cl')~a{Ci)-a{R^) 

= \UI -foO f^Ui)] + \ipQ o ip\ui) - (Po{Ur)\ - \ui - (po{Um)\ - |¥'o("r) " ¥'o(«m)l 

= -2|^^(«,) - <p^oium)\ < -2Lip(y,^)a(i?i). 
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Case CC-1 : (C|C^)-(C^'). This is another simple case. We begin with 

max <p*(um)) <Ur < Um < Uh and < Um- 

When < Ur, all the waves are classical and it is easy to show that [V] — 0. When 
Ur < 0, we still have >fi^{um) < Ur and therefore the important states are ordered 

(foiUr) <Um<Ul. 

It then easy to check that even in this case, [V] = 0. 

Case CC-2 : {C^C^)-{C^'). This interaction is constrained by the initial states in 
the following manner 

(^'(U;) <Um <Ur < (^"(Wm) < Ul, and Um < 0. 

Two subcases appear depending on the sign of Ur- When Ur > 0, the interaction is 

{C\.C^)-{C^ ). Since (p^{ui) < Um, '^oi'^m) < ° f^i'^i) < Combining this with 
the fact that Ur < (p^{um) < ¥'o(^to)> can deduce that the states appearing in [V] 
are ordered 

Ur < iPoiUm) < Ul. 

The change in V is then 

[V]=a{Ci)-a{Ci)-a{Cl) 

= \UI - Ur\ - \UI - ipQ{Um)\ - \(pQ{Um) - Ur\ =0. 

The second subcase treats Ur < and interactions (C^Cl)-(Cj^ ). The states used 
in our definition of the strength of the waves are 

(foiUr) < iPoiUm) < Mi- 

A short computation shows that 

= \UI - ^oM\ - \UI - Vo("m)l - \V>o{Um) " <^o('"r)| = 0. 

Case CC-3 : {C^C^)-{Nj_ C^' ). This interaction represents a typical transition 
from one crossing shock to a nonclassical shock. The states are 

(p^{ui) < (p\um) <Ur < <P^{ui) < Um < Ul, < Um- 

The first (and most common) subcase occurs when Ur < and the interaction is 
{C\_Cj.)-{Nj_ Cl_ ). The first subcase needs to further subdivided into two cases. 
When Um < fo ° 'P^{ui), all the important states arc ordered 

(PoiUr) <Um<<fo° f^'^l) < ^l- (3.14) 
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Under these circumstances, the functional V doesn't change. 
[V] = a{Ni) + a(Cl') - a{Ci) - a{Ci) 

= \ui -(PqO (p\ui)\ + \lPq o ip\ui) - (po{Ur)\ - \ui - Um\ - \Um - ¥'o("r)| = 0. 

On the other hand, when ip^ o ip^{ui) < Um in the first subcase, we have 

Lp\^{Ur) < Lp\i O 1^^ [Ul) < Um < Ul. (3.15) 

It is easy to see that we again have \y] =0. 

In the second subcase, < u^, we again need to introduce two additional subcases to 
handle the interactions {C\_C\.)-{N^ ). When < tp^ o (p^{ui), then the states 
used in the definition of wave strengths are 

Ur <Um < VoO V^iui) < Ui- 

The change in [V] is 

[V]=a{Ni)+a{Cl)-a{Ci)-a{Ci) 

= \UI - ° <^^"/)l + ° ^Hui) - Ur\ - \UI - Um\ - \Um - Ur\ = 0. 

When < Ur and tpg o ip^{ui) < Um, we obtain the same result. 

Case CN-1 : (C|Afi)-(Ci ). The states defining the waves are characterized by 
the inequalities 

0<Um<ui and «r = ¥'''(Mm). 

This implies that 

</?0 ° 'P\um) <Um< Ui. 

We deduce 

[V] = a{C^^)-a{Ci)-a{Ni) 
Case CN-2 : (CiA^^)-(C|'). We begin with states satisfying 

ip^{ui)<Um<0 and Ur = ip'{Um)- 

To measure the wave strengths, we observe that ip^{ui) < Um implies 
<Po(^m) < ¥'o ° < ^iid therefore 

The jump in V is now 

[V]=a{Ci)-a{Ci)-a{Nl) 

= \UI - <f\um)\ - \UI - >fo{Um)\ - \fo{Um) " f\um)\ = 0. 
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Case CN-3 : {C^N^)-{Ni' Cl'). The states are initially ordered as 

< Um < ui and Ur = tp^{Um) < ^^Ui). 

A first subcase occurs when Um < Vo ° v'^iui). To compute the change in V, we can 
then use the inequalities 

fO ° f^Um) <Um<<PQO ip\ui) < Ul, (3.16) 

to deduce 

[V]=a{Ni)+a{Cl_)-a{Ci)-a{Ni) 

= \UI - </?o ° + \Vq ° f^Ul) -(PqO (p\um)\ 

- \ui - Um\ - \Um -^0° {Um)\ 
= 0. 

Similarly, if ip^ o ip^{ui) < Um, then the important inequalities become 

<fiO ° V>\um) <^o° V^Ul) <Um< Ul, (3.17) 

and [V] = 0. 

Case NC : {N^C'' )-{C^'). This interaction is constrained by the states 

Um = 'P^{ui) and (p'^{ui) <Ur < V5*(Um) < Ul. 

The first subcase occurs when Ur < and the interaction is precisely {N^C^_)-{Cj_ ). 
The important states are then ordered as 

foiUr) < (PoO(p\ui) < Ul. 

With these observations, we find that 

[V] = a{Ci')-a{Ni)-a{Cl) 

= \UI - </?o("r-)| - \UI - <^o ° V'^wOI - l<^0 ° V>\ui) - (Po{Ur)\ = 0. 

In the second subcase, Ur > 0, it is easy to check that 

Ur < tp^O (p^{ui) < Ul, 

and [V] = 0. 

Case NN : {N^N^)-{C^ ). This interaction is the limiting case Ur — > f^ium) of 
Case NC. By continuity of wave strengths, we must also have [V] = 0. 

□ 

4. Quadratic interaction potential (part 1). In the rest of this paper, we 
investigate two quadratic interaction potentials, keeping in mind from experience 
with classical shock waves, that different functionals may be of particular interest in 
different circumstances. We begin by searching for a functional of the form 

Q{u{t)):= Yl cr{ur,u^)a{u?,u?), (4.1) 

a approaches 
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where the proper definition of "pairs of approaching waves" is essential and is now 
specified. 

In GHmm's original paper for systems of conservation laws, a definition was 
proposed which, in the scalar case, amounted to stating that pairs of waves were 
always approaching unless both are rarefactions. The purpose of this section is to 
investigate the original definition of Glimm in the context of nonclassical shocks. 

Definition 4.1. A wave a is said to weakly approach a wave B, unless both are 
rarefaction waves. As far as this definition is concerned, waves R\_,R_ are both to be 

I I T T T 

considered as rarefaction waves are all to be considered 

as shock waves. 

Our main result in the present section is as follows. 

Theoreim 4.2 ("Weak interaction" potential for nonclassical shocks). Let ip^ be 
a kinetic function satisfying the properties (A1)~(A3). Consider the functional Qweak 
defined by (|4.ip where the summation is made over all weakly interacting waves in 
the sense of Definition \4-l\ Then, when evaluated on a sequence of front-tracking 
solutions, Qweak is strictly decreasing during all interactions except in the cases RC- 
3, CR-4, CC-3 and CN-3. In fact, for each of these exceptional interactions, there 
exist initial data for which + Cq Qweak is increasing during the interaction for every 
positive Cq. 

In contrast, in Chapter 8 of [S] in joint work by Baiti, LeFloch, and Piccoli, 
different definitions of both wave strengths and approaching waves are used and the 
resulting Glimm functional V+K Q is strictly decreasing for some K. In this sense, the 
interaction functional Qwoak above may appear to be less satisfactory. However, our 
assumptions on the kinetic function are completely natural -a major advantage toward 
a future extension to systems- and, furthermore, an analysis of "splitting/merging" 
solutions (in the following section) will show that globally in time the functional Qweak 
does decrease. 

Several justifications for our definition of potential are now provided, the strongest 
argument being the requirement of continuity: 

1. Given that V is continuous in BV(M) and that 'I'q is already an isometry with 
respect to this BV(IR.) norm, it is tempting to assume that any reasonable 
interaction potential Q should also be continuous in BV(R). We observe that 
any shock C\_ can be continuously deformed (as measured by Definition 12. ip 
into, first, a crossing shock Cj_ and, then, a pair of shocks N^Cl. It IS easy 
to see that imposing continuity would imply Definition 14.11 

2. Another argument can be found by looking at a class of solutions called split- 
ting/merging solutions, introduced in [9] and discussed further in Section [5l 
These solutions illustrate that some initial data can go through a nearly peri- 
odic process of creation and destruction of nonclassical shocks. In particular, 
nonclassical shocks can indirectly have non-trivial interactions with shocks 
on their right-hand side. 

3. In [6j, nonclassical shocks are precluded from interacting with their right-hand 
neighbours, and it is argued that nonclassical shocks are (slow) undercom- 
pressive and, thus, move away from their right-hand neighbors. However, 
this definition of approaching waves ignores the above-mentioned possibility 
of nonclassical shocks having indirect interactions with shocks on their right- 
hand side. In any case, such an interaction functional then would not be 
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continuous in BV(R). 

For the proof of Theorem 14.21 we proceed as follows. During any isolated inter- 
action between two waves in an approximation, the change in Q — Qwcak is of two 
types 

[Q] = [Q]i + [Qh- (4.2) 

In this decomposition, denotes the change in the products of the strengths of 
waves either incoming or exiting the interaction and [Q\2 denotes the change in prod- 
ucts of strengths of waves where only one of the waves was directly involved in the 
interaction. Moreover, if a wave C is involved in an interaction, we define 

W{C) := 

B approaches C 
B did not interact 

According to Definition l4.1[ if the incoming wave Cin and the outgoing wave Cout are 
of the same type (i.e. both of rarefaction or shock type), then W{Cin) — VF(Cout)- 
Theorem 14.21 requires the following lemma which emphasize the (wrong) sign of the 
terms [Q]i appearing in the most difficult interactions. The lemma gives a precise 
characterization of the failings of Qweak for those interactions. 

Lemma 4.3. The interactions involving a crossing shock, classical or not, and 
a small wave Win are RC-3, RN, CR-4, CC-3 and CN-3. Suppose that the states 
neighboring the two consecutive waves are ui , Um and Ur ■ For ui fixed and as a function 
of the strength of the incoming small wave, the change in [Q] i is as follows : 
RN : [Q]i < for all 

CN-3 : [Q]i < if Cin is sufficiently weak; 
RC-3 : [Q]i < if Rin is sufficiently strong; 

CR-4, CC-3 : [Q]i can take on both signs, depending on Um and Ur- 
Moreover, for RC-3 and CR-4 interactions, the largest positive change occurs when 
a{R,n) = 0. 

Proof. We distinguish between several cases. 
Case RN : {r\n^)-{N^ R[_ ). The states appearing in the shock strengths 
describe two subcases: 

either i^q ° V^i^i) < ^p' {um) < ui < u™, 

or (fl O (p^{ui) < Ul < IPqO Lp'{Um) < Um- 

For both sets of inequalities a{R[_ ) < cr{R\_) and cr{N^ ) < cr{N^). As a result, 

[Q]i = a{Ni')a{Ri') - a{Rl)a{Ni) < 0. 

Case RC-3 : (i?^C^)-(A^j_ C^'). The fundamental states describe two subcases 
(13. 7p and p.8p . here rewritten 

either ipQ{ur) < V'o ° V^M < < "m, 

or Ur < IfQ O ip^{ui) < Ui < Urn- 
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Focusing on the first case, the change in [Q]i is 
[Q]i{ui,um,ur) =a{Ni')a{C^') - a(i?|)a(Ci) 

= (u; -(foO (p^{ui)) {(fiQ O (p\ui) - (po{Ur)) - {Um - Ui) {um - 'Po{Ur)) ■ 

Clearly, the RC-3 interaction can still occur even as the strength (T(i?^) 0, 
therefore the expression above shows that [(5]i > in that limiting case. In fact, we 
will show that this is the largest value of [Q]i- We begin by computing 

= - {Um - V>oiUr)) - {Um - Ui) < 0. 

Then we evaluate [(5]i when Um = p{ui), where 

p{ui) := ui + (ui -iPqO (p\ui)) =ui + a{N^'), 

and find 

[Q]l{ui,p{ui),Ur) ={ui - ipQOip^Ul)) {(fioO(fi\ui) - <fio{Ur)) 

- {ui -ifioo <fi\ui)) {2ui - ° <p''(wi) - VoiUr)) 
= - 2{ui - o <p\ui))^ = -2a{Ni'f. 

Therefore, for the unbounded set of states satisfying Um > p{ui) (i.e. 

a{Rl) > (T(iV±')), we have [Q]i < 0. 

For fixed ui , the remaining set of states describing an RC-3 interaction form a 
bounded set in the {um, Ur) plane defined by 

Um < P{ui), 

Along the edge Um = p{ui), [Q]i is negative and since [Q]i is decreasing with respect 
to Um, the largest value must occur along the edge Ur = ip^{u„i)- We therefore 
compute the change in [Q]i along that edge, using u^ as it's parameter. 



-f— ([Q]l{ui,Um, ^p\um)) 
(iUui V 



- (Um - Ul) {Um - </Jo ° ^Hum)) ^ 
= in-^^,0^\uO){-^^M) 

- (u„j - (foO tp'^{Um)) - (Um - Ul) ^1 - '^^^^^^ {UrrSj 

Using the inequalities for the states and the fact that 

Lip((^Q o (^f) < Lip((pg o ip^) < 1, we deduce that the terms above are negative and 
therefore, that [Q]i must attain it's largest value for the smallest value of Um, 
namely when Um = ui = f^iur) and (t{R\_) = 0. 
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Assume now that (p^{ui) > and that Ur could be positive but continue to fix ui. 
Then the set of states describing an RC-3 interaction form an infinite strip in the 
{um,Ur) plane defined by 

Ul < Um, and V^{ui) < Ur < ^'^{ui), 

where we have assumed the reasonable fact that {(p^)' > at = f'^{ui) (in fact, 
since {(p"^)' < at the origin, we want to assume only one change of sign for (lys")' 
and more general cases are similar). The same calculation as above shows that 
d[Q]i/dujn < and therefore that the maximum value of [Q]i must occur when 
w„i = Ul, i.e. cr(i?i„) = 0. 

Case CR-4 : {CiRi)-{Ni'cl'). The states satisfy (j?!^ 

'PoiUm) < VoiUr) < Ifi^ O P^ M < Ul, 

and therefore 

[Q]iiui,u^,u,) =a{Ni')a{Cl_') - a{Ci)aiR^) 

= {ui -ipoo (p\ui)){'p^a ° - VoM) 

- {ui - p\i{Ura)) (Vo(^r) - </5o(^m)) ■ 

If we let Ur and Um approach each other while maintaining the condition 

Ur < ^'^{Ui) < Um, 

then ip\^{ur) — (Pq(uto) — > while ui — (pQ{um) remains bounded. Upon inspection, it 

I / I / 

is clear that in this limit, [Q]i ). We now demonstrate that this is 

the largest value of [Q]i- When ui is fixed, the admissible states belong to a 
rectangular domain 

{Um,Ur) e [p'^iui), 0] X [p\ui),ipi{ui)]. 
With respect to Um and Ur, we find 



d[Q]i _ dip'o 
dum du 



(Urn) (^{ipliUr) - ipliUm)) + ("i - V3o(^m))) < 0' 



and 



d[Q]i _ dip', 
dur du 



(Ur) (^{ui - PqO ip^{ui)) + {ui - (/?o(Um))) > 0. 



Taken together, these inequalities imply that the largest value of [Q]i occurs at the 
corner of the domain where Um — Ur — p^^iui), i.e. when a{R^_) — 0. On the other 
hand, the smallest value would occur along the boundary Ur = (p'{ui), where one 
finds a negative quantity 

[Q]i ={ui -VqO ip\ui)){(pl o (p\ui) - (Po o 'P^Ul)) 

- {ui - p^,{u„i)){iPo O p\ui) - PoiUm)). 
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Case CC-3 : {C\^C-^)-{Ni' C^'). As mentioned in the proof of Theorem O the set 
of states are subdivided into two subcases depending on the sign of Ur- When 
Ur < 0, we identify two possibihties (|3.14p and (|3.15p . namely 

either </?o(wr) < Um < {ui) < uii 

or (p\i{Ur) < V'O ° V^i^l) < Urn < Ul, 

where the second set of inequahties corresponds to a weak C\_ . The change in [Q] i is 
then 

[Q]i{ui,u^,ur) =a{Ni')<j{Cl') - <j{Ci)<7{Ci) 

- {Ul - U„i)(u,„ - ipQ{Ur)). 

Consider the function B{X) :— A(l — A) whose maximum on the interval [0, 1] is 
attained at A = 1/2. There exists constants C, A and A' such that the change can be 
rewritten as 

It is therefore clear that will be negative if and if |A — 1/2 1 < |A' — 1/2 1. 
Unfortunately, this set cannot be described in a simple manner in terms of the 
strengths of the waves. 

Case CN-3 : (C| A^i)-(iVi'cl'). We identify two subcases 

either ip^ o ip^{um) < Um < (Po o p^{ui) < ui, 

or Ifl O ip^{Um) < iPo° <Um < Ul, 

where the second one occurs when the incoming shock C:{i is weak. The change can 
be written as 

Using a bit of algebra, we rewrite the change as 

[Q]l ={ui - <^o ° V^M){^0 ° - ipoO ip^{Um)) 

= - {um - V'o O V^iui)) {{Ul - Um) - (v'o ° M " Vo ° (4-3) 

thus concluding that [Q]i < as long as C\_ is weak. 

□ 

Proof of Theorem \4-^ Throughout the proof, we use liberally the estimates 
derived in the proof of Theorem 13.11 
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RC-1 : iRlC^)-{C^'). After examining p.Sp . it is immediate that 
cr(C^) > cr(C^') and W{C^) = V7(C^'). Therefore, 

[Q]i = -a{Rl)a{C^) < 

and 

[Q]2 = W{C^')a{C^') - W{Rl)a{Rl) - M^(C^)ct(C^) 
= -WiRl)a{Rl) - W{C^){a{C^) - a{C^')) < 0. 

Case RC-2 : (i?^Cj_)-(A^j_ R{_ ). The conditions defining this case imply two 
subcases : 

either ip^ o (p^{ui) < ipo{ur) < ui < Um, 

or Ifil O (p^{ui) < Ui < tPoiUr) < Um- 

Given that ip^Um) < Ur < 0, then property (|2.ip leads to 
This means that in the first subcase, 

cr(i?i ) = l-^ol^r) - "^0 ° < Ifo ° f\u„i) - iPqO (p\ui)\ <\um~Ui\^ 

In that subcase, we also have 

a{Ci) ^ aiRl) + \ui - ifiliur)], 
a{Ni')=a{Ri') + \ui-ip\,{ur)l 

then cr(-/V± ) < (j(CJ^). Similar arguments prove the same inequalities in the second 
case. 

These inequalities therefore imply that 

[Q]i = a{Ni')<j{R{_') - a{Ci)a{Rl) < 0. 
Using the proposed definition of weakly approaching waves, we also deduce 

[g]2 =W(Ni')(j{Ni') + W{R[_')cr{R[_') - W{Rl)cr{Rl) - W{Ci)a{Ci) 
= - W{R{_'){(j(rI) - (t{r[_')) - W{Ni') {(j{Ci) - (j{Ni')) < 0. 

Case RC-3 : {R^_^C^)-{N^ C^'). Two subcases occur depending on the sign of Ur- 
If Ur < 0, then (|3.7[) holds and when Ur > 0, then (|3.8[) holds. In both of these cases 

aiC^)~cT{Rl)^cT{C^') + a{Ni'). (4.4) 

In Lemma 231 we already showed that [Q]i could be positive. For the other 
interaction term, we use (|4.4p to verify 

[Qh ^W{Ni')aiNi') + WiC^'MC^') - WiRl)aiRl) - V7(C^)a(C^) 
- - W{Ni') {a{C^) - a{Ni') - <j{C^')) - W{Rl)a{Rl) < 0. 
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Case RN : {RlNi)-{Ni' rI'). The change [Q]i was studied m Lcmma|131 For 
the second term, using the bounds on the wave strengths given in the previous 
lemma, we again have a negative contribution 

[Q]2 =WiNi')a{Nl') + W{R[_')a{R^') - W{Rl)a{Rl) - W{Ni)a{Ni) 
= - W{Ni) {a{Ni) - cr(iVi')) - W{rI) {a{Rl) - cr(i?i')) < 0. 

Case CR-1 : (C^i?i)-(Cj_ ). In this case, the states satisfy 

foium) < foiur) < 'Po{fH''^i)) < ^he wave strengths satisfy (t(C^ ) < o-(Cj_). 

Since only one wave is outgoing, [Q]i = — cr(CJ^)cr(i?i) < 0. On the other hand, it is 
easy to check that 

[Q]2 ^W{Ci')a{Ci') - W(Ci)a{Ci) - W{R{)a{R{) 
= - W{Ci) {(j(Ci) - (j{Ci')) - W{Ri)a{R^) < 0. 

Case CR-2 : {C\_R^_^)-{C^ ). This case is entirely classical so it is easy to check 
that 

[Qh^-<j{Ci)a{Rl)<0, 

[Qh = - WiCi) {a{Ci) - aiCi')) ~ < 0. 

Case CR-3 : (Cii?i)-(7V4'i?i'). A first subcase is defined by the additional 
condition ui < ip^{ur) which provides p.lOp 

<Po("™) < Vq° V^{ui) <ui < (filiur). (4.5) 

The relative strengths of the waves are then seen to be (t{R^_ ) < ^(i?!) and 

<t(-/V^ ) < a{C^). In the second subcase, given by p. lip , these two inequalities still 
hold. It is now easy to conclude 

[Q]i ^a{Ni')a{R{_') - a{Ci)a{Ri) < 0, 

[Qh = - W{Ci) {a{Ci) - a{Ni')) - W{R^) {a{R^_) - f7(i?i')) < 0. 

Case CR-4 : {CiR^)-{Nl' Cl'). The inequahties (|3T3| defining this case suffice 

to show that cr(A'^i') + cr(Cl') < a{Ci). In Lemma l43l we showed that [Q]i > 
but the second term is negative 

[Qh = -W{Ci) {a{Ci) - a{Ni') - a{Cl')) ~ W^(i?i)a(i?i) < 0. 

Case CC-1 : {C]_C^)-{C^'). When Ur is positive, all waves are classical and it is 
easy to show that [Q]i < and [(5]2 = 0. It is an exercise to see that when Ur <0, 
then [Q]i < is still negative and [Qh vanishes. 

Case CC-2 : (c4cT)-(Ci'). Two subcases appear depending on the sign of Ur but 
each time, 

[V] = a{C^') - a{Ci) - a{C^) = 0. 
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Only one wave is outgoing, so [Q]i < and the previous identity implies that 

[Qh - 0. 

Case CC-3 : {C\.C^)-{nI' C^' ). Since all waves involved are shocks and [V] ^ 0, 
then [Q]2 = 0. 

Case CN-1 : (C|7V4)-(CJ^ ). This is another simple case where the fact that 
[V] ~ and that W{-) is equal for all waves involved, implies [Q]2 = 0. 

Case CN-2 : (C^ A^^)-(C|'). Same as Case CN-1. 

Case CN-3 : (C|iVi)-(Afi'cl'). Two subcases occur but each time [V] = 0. All 
waves are shocks so [V] = implies that [Q]2 = 0. 

Case NC : {NiC^)-{C^'). Same as Case CN-1. 
Case NN : (iVi7V]:)-(C|'). Same as Case CN-1. 

□ 

5. Global diminishing property for splitting-merging patterns. In this 
section we show that, despite the fact that the quadratic interaction potential Qweak 
increases during interactions of type CR-3, RC-4, CN-3 and CC-3, this potential 
is indeed strictly decreasing globally in time for a large class of perturbations of 
crossing shocks. Hence, with the local bound we describe below, we provide the 
first steps towards an analysis of the global-in-time change of Qwoak for arbitrary 
nonclassical entropy solutions. This section, therefore, provides a strong justification 
for the potential proposed in the previous section. 

The splitting-merging solutions considered now were introduced in LeFloch and 
Shearer [9], where a modification of the total variation functional [Bj was shown to 
be strictly decreasing along the evolution of such splitting-merging solutions. The 
total variation functional V presented in Section [3] also accomplishes this, but here 
we improve on those results by establishing a similar monotonicity result for the 
quadratic functional Qwcak- Our analysis also brings to light some interesting aspects 
of splitting- merging solutions that were not seen in [9] . 

Splitting- merging solutions are, roughly speaking, perturbations of crossing shocks 
that lead to the creation and destruction of a nonclassical shock. Such solutions con- 
tain two (classical and nonclassical) big waves that may merge together (as a clas- 
sical shock) and also interact with (classical) small waves. A typical initial data for 
splitting-merging patterns is formed of 

(i) an isolated crossing shock with left- and right-hand states satisfying 
w_ > and (p^U-) < u+, but ii+ — ip^u-) small, 

(ii) followed, on the right-hand side, by a small rarefaction and a small shock. 
The rarefaction is sufficiently strong that it has an interaction of type CR-4 with the 
crossing shock, thereby leading to the creation of a pair of shock waves Ni,CL If the 
right-most shock is sufficiently strong, then when it eventually interacts with and 
the resulting shock will begin to approach N\_ . The final interaction of type NC will 
involve Nj_ and the shock just described, thereby eliminating the nonclassical N^. By 
adding more waves to the left and the right, this process of creation and destruction 
of can be repeated indefinitely. 

We consider a slightly more general configuration in the sense that we do not 
explicitly demand that a small shock on the right be responsible for the penultimate 
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NC interaction. Fix some value u* > and define 

Unix) = < u, , (5.1) 

°^ ' [ipHu*), x>0. ^ ' 

Let be some function of locally bounded total variation and of oscillation bounded 
by some small positive e, i.e. 

a{0,{x),O) < e, X eR. 



Without loss of generality, we may assume that 9e is piecewise constant. Let u* be 
the nonclassical solution to the conservation law (II. 1|) with initial data Uq + 9e, as 
generated by the front-tracking method. Assuming the solution initially possesses 
a single isolated crossing shock located at the origin x — 0, that is, assuming that 
w*(0— ) > 0, we see that the crossing shock will be adjacent to many small classical 
shocks and rarefactions. After an interaction of type RC-3, CR-4, or CC-3, the small 
waves neighboring may lead to the creation of a pair of waves, a nonclassical shock 
Nj_ and a classical shock . After the creation of Nj_ , the only types of interaction 
involving small waves incoming from the left of are RN and CN-3. The only types 
of interaction involving small waves and the shock , coming from either the left or 
the right, are entirely classical (CC-1, RC-1 or CR-1). Moreover, no waves can cross 
from the right or the left although the small waves that reach Nj_ from the left, will 
cross and eventually reach . Therefore, the only way that the nonclassical shock 
Nj_ can be destroyed is if the shocks Nj_ and change their speeds and eventually 
interact back together, leading us back to (a perturbation of) the original crossing 
shock. 

These observations, in particular, imply that no waves can exit the domain Q 
bounded by the trajectories of A''^ and . Our goal, in the present section, is to 
obtain a local bound on the change of the potential Qwcak relative only to the waves 
entering the domain Q. It should already be clear that the key here is comparing the 
total strength of the waves crossing to the total strength of the waves terminating 
at CT. 

Before stating our main result, we introduce some further notation. Let to be a 
time of creation of a nonclassical wave N^. and denote by ^1,^2, • ■ • ,t,n the times of 
the next m interactions between Nj_ and small waves Wi on the right, and let t„i+i be 
the time at which Nj_ is destroyed from an interaction with the shock . Similarly, 
let ti and ti be the times at which an interaction occurs between and the left 
incoming waves Wi or the right incoming waves Wi , respectively. We define the total 
variation along the trajectory Nj_ to be 

m 
i=l 

and its signed variation 

SY{Ni) := a(A^i(Wi-)) - a{Ni{to+)). (5.3) 
Completely similar definitions also apply to the wave , but an additional decompo- 



(5.2) 
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sition exists by separating the contributions from the left- and the right-hand sides: 

= :SVl(CT)-^SVk(CT), 

where Si (s^) is +1 if Wi (Wi) is a shock and —1 otherwise. For convenience, the 
strengths of the smaU wave Wi, before and after it has crossed Nj_ at some time U, 
are denoted by W~ and , respectively. 

Theorem 5.1 (Global diminishing property for splitting-merging patterns). Let 
ip^ be a kinetic function satisfying the properties (Al)-(A3), together with the mild 
requirement at some given reference state m* > 0. 

L* := {ip'^o^')'{u*)e {1/2,1). 

Let u* he the nonclassical solution to the conservation law (II. 1[) with initial data 
Uq -\- 9^, where Uq is defined in (15. 1|) and the perturbation 0e is of locally bounded 
variation and of small amplitude, that is, ||0e||L°°(R) < £■ Suppose that u* exhibits a 
splitting-merging pattern on the time interval [to;^m+i] along successive interactions 
with small waves (cf. the notation above). Provided e is sufficiently small and the 
total effect of all waves on the classical shock increases its total strength, that is, 

SV(CT) > 0, 

then the variation of the potential Q — Qweak within the domain Q bounded by the 
trajectories of Nj_ and is negative 

[Q]iL<0- 



Recall that the expression [Q]i was defined earlier in Section 4. In the above 
theorem, only contributions within the domain are considered. Two important 
remarks about our assumptions are in order: 

1. The condition on L* is satisfied for a wide range of fluxes, entropies and 
nonlinear diffusion-dispersion models, see Theorems 2.3 and 3.3 in [B]. In 
particular, it is true for the cubic flux f{u) = u^—au, the entropy U (u) = u^/2 
and all ratios of diffusion over dispersion. The non-trivial aspect of this 
condition is the lower bound L* > 1/2 since the upper bound is always 
satisfied. 

2. Our main assumption SV(C^) > requires that the total effect of the inter- 
action of all waves on is to increase its strength. In fact, this is always the 
case for the perturbations of splitting-merging solutions within the setting 
[S]. In our slightly more general setting though, waves crossing through Nj_ 
will change the critical state ip^{ui) and therefore could conceivably lead to 
an NC interaction even if interacted only with rarefactions on the right 
(SV(C^) < 0). Intuitively, one would like to show that the existence of an NC 
interaction at time tm+i implies that the wave C became strong enough to 
change course and therefore that SV(C^) > is satisfied. This is the subject 
of ongoing research. 
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It is interesting to check that the condition SV(C^) > of Theorem lS.ll is in fact 
necessary, by studying a simple situation involving a single rarefaction R crossing the 
nonclassical shock Nj_, and one classical shock C interacting from the right with 
(assuming that (^^ is monotonically decreasing, say). Suppose that the strength of 
the rarefaction wave after crossing is i^r and that the change in during this 
interaction is vl. In Lemma 1 5 . 2 1 b elow . we will check that ly^ < i^r and that + t^R 
is the strength of the original rarefaction. So, assuming that C" interacts with 
before the rarefaction does, then the total change in Q is given by : 

= +<j{Ni)(7iC^) (duetoCR-4), 

- iyL{(T{N^) + i^b) (duetoRN), 
-(j{C^)a{C') (duetoCC-1), 

- (cr(CT) + a{C'))i^R (due to RC-1), 

- (a(iVi) - i^l) (a(CT) + a(C') - i/fl) (due to NC). 
The last term can be rewritten as 

-<j{NiMC^) - a{Ni){a{C') - ,^r) + z^L(a(CT) + a(C')) - i^LiyR. 

The contribution from the RC-1 interaction can control the third term. Unfortunately, 
if fR is large with respect to cr(C"), i.e. SV(C") < 0, then the CC-1 term cannot be 
used to control the second term. 

For the proof of Theorem lS.ll we will need the following two lemmas. The interest 
of the first lemma is to make more precise the (mainly linear) dependence of the 
change [Q]i in terms of the incoming wave a{W~). The second lemma provides an 
estimate which closely relates the signed variations of and iV^ . 

Leimivia 5.2 (Interactions with the nonclassical shock). Consider interactions RN 
and CN-3 at the time ti involving a left-incoming weak wave Wi and the nonclassical 

4- 



shock Nj_ . Then, there exists a positive constant Li < 1 such that 



a{W+) = LMWrh 
a{Niit,+)) = a{Ni{t,-)) + s,(l - L,)(t(T4^-), 

[Q]i = -(1 - L,)a{W-) (a(7Vi(t,-)) - s.LMWr)'^ < 0, 

where Si is +1 if Wi is a shock and —1 otherwise. Moreover, one has 

L, = L* + 0{e). (5.4) 



Proof. We consider only a RC-3 interaction, since the calculations for an incoming 
rarefaction are similar and have been essentially treated in Lemma l4.3[ equation (j4.3[) . 
When Wi is small, then the states are ordered as in ()3.17|) . namely 



Then, we have 



< u 



a{W^^ 



L,\u] 
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with Li < 1 since (p^ oip^ is a contraction. On the other hand, relation (15. 4p is obvious 
since the u\,i,u\ belong to a neighborhood of u* of size e. 
Finally, the outgoing nonclassical shock has strength 

while the change [Q]i takes the form 

[Q]i = a{Ni{U^))a{W+) - a(7Vi(i,-))a(W,-) 

= (cy{Ni(U-)) + (1 - L,)^(W,-)) i.a(Tyr) - ^{Ni{U-))cy{W-) 
= -d - Li) ^{W-) (cj{Ni{U-)) - L,cj{W-)) < 0. 

□ 

Lemma 5.3 (Property of the signed variations). Setting A* := (1 — L*)/L* , one 

has 

< 0(e)TV(7Vi). (5.5) 
Proof. Using Lemma 15.21 we compute 

rn 

SV(7Vi)=^s,(l-L,)'T(W,-) 

m rn 

= J2s,{l- L*)<j{Wr) + ^ s,iL* - LMWf) 

711 rn 

=X* s.L*a{Wr ) + ^ s,{L* - L,;)a(W,-), 

thus 

rn 

SV{Ni) = A* SVl(CT) + (1 + A*) ^ s,(L* - L,)a{W-). 

1=1 

For each index i, L* — Li = Oi{e), so if we take 0(e) = maxiOi{e) and use L := 
mini Li, then the difference is 

m 

<20(e)5^|a(W,-)| 

i=l 

< yO{e) TV(iv4) = 0(e) TV(7Vi). 

□ 



(^^o^^(«J)-^^o^^«) 



A*SVl(CT) -SV(iVi) 



A* SVl(CT) -SV(iVi) 
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Proof of Theorem \5.1l To help the reader understand the wave interactions (and 
cancellations) in this proof, we begin with a few preliminary comments. Simply put, 
the final NC interaction should provide a quadratic term which is, later, cancelled by 
a similar quadratic term when the two waves merge. One may naively expect the sum 

-a{Ni{t^+i-)) a{CHt„,+i-)) + cT{Ni{to+)) a{CHto+)) 

to be negative. Of course, the cumulative strength of the changes during the in- 
teractions with the small waves must also be taken into account. The difference 
between the strength of the initial and final waves Nj_ , is measured by the signed 
variation SV along those two shocks. Our proof below shows that, along the trajec- 
tories and , the change [Q] i is negative and proportional to the total variation 
TV(iV^) + TV(C^). The total variation being larger than the signed variation, after 
further analysis one can conclude that [QJijn < 0. 

The key technical information is provided by Lemma 15.31 which implies that, 
up to a quantity of order ©(e) TV(iVi), we have SV(7Vi) = A* SVl(CT) and, in 
particular, that the signed variations have the same sign. We now have all the tools 
necessary to proceed with the proof of Theorem 15.11 

The perturbation 0^ has bounded oscillation and therefore it can only alter the 
right-hand state of by an amount e. The small waves entering fl through only 
alter its strength by |(1 — Li) a{W~)\ < e. In both cases, we expect that everywhere 
along their trajectories, 

a{Niit))>a{Ni{to+))~2e, 
a{CHt)) >a{CHto+)) -2e. 

Ignoring the negative contribution to [Q]i coming from the interaction that gen- 
erated Nj_ (which, anyway, can be arbitrarily small), and neglecting also all classical 
interactions inside H. (for which [Q]i < and possibly 0), we have 

[Q]i\^< + <7{Ni{to+))<j{C^{to+)) 

rn 

- ^(1 - LMW-) {<j{Ni{U~)) - s.LMWr)) 



i=l i=l 

-a(7Vi(Wi-))a(CT(t„+i-)), 
and after simplification 

[Q]i[^< + <j{Ni{to+))a{CHto+)) 

- (a{Niito+)) - 3e) TV(iVi) - (a{CHto+)) - 2e) TY{C^) 

- ((T(iVi(to+)) + SV(iVi)) (a(CT(io+)) + SV(CT) 



Introduce constants tjn, rjc & [0, 1] to be determined later in the proof. The terms 
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can then be split in the foUowing way. 

<-((!- VN)a{Ni{to+)) - 3e) TV(7Vi) (5.6) 
- ((1 - r^cy{C^to+)) - 2e) TV{C^) 
-a{Ni{to+)) (?7jvTV(7V4) + SV(CT)) 
-a{C^to+)) (r,cTV(CT) + SV(iVi)) -SV(iVi) SV{C^). 

If riN, Tjc are known a priori, the first two terms are negative (and can be neglected) by 
taking e sufficiently small. Clearly, the main difficulty now lies in the sign of SV(iV£) 
and SV(C^). If both are positive, then all the terms are negative and we have proved 
our result. Recall that SV(C^) > 0, by assumption, so that it suffices to consider the 
(only possibly unfavorable) case 

Lemma [5.3l implies that we can assume that SVl(C^) is also negative, provided we 
accept a small correction term 0{e) TV(iV^). From this, we deduce that SVfl(C^) > 
and in fact, necessarily 

|SVl(CT)| < |SVj^(CT)|. 

We immediately note that any correction terms of order 0{e) TV(iV^) can be included 
into the first term of the decomposition (|5.6[) and, therefore, taking a smaller e if 
necessary, we can make the new term negative: 

-((1 - TlN)(y{Ni{t^+)) - 0(e)) TV(iVi) < 0. 

This fact is used below without further comment. 

Using Lemma [531 we have |SV(A^i)| « A*|SVl(CT)| < A*TV(CT), and the 
fourth term in (15.61) can be written as 



-a{CHto+))(TicT\{C^) + SY{Ni)^ 

<-a(CT(to+))(r/cTV(CT)-A* TV(CT)) + 0(e) TV(7Vi). 

Taking rjc = A* suffices to guarantee that the term vanishes. We now re-organize 
the third and fifth terms of (|5.6p to take advantage of the known signs of the signed 
variations, as follows: 

A:=- a{Ni{to+)) (tjn TV(iVi) + SV(CT)) - SV(A^i) SY{C^) 



<y{Ni{to+)) (rjN TV(iVi) + ^ SVl(CT) 
- <y{Ni(to+)) ((1 - ^) SVl{C^) + SVr{C^) 



Vn 

2 ' 

- SV(iVi) SVl(CT) - SV{Ni) SVfl(CT) + 0{e) TV{Ni). 
Observe that we have the following lower bounds 

(1 - ^) SVi(CT) + SYniC^) > ^| SYUC^) 
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and upper bounds 

SVl(CT) =^ SV(iVi) +0{e)TV{Ni) <2TV{Ni) + 0{e)TY{N. 

A 

Using the two previous bounds, we find 

A<- a{Ni{t„+)) TV(iVi) - m TV(iVi)) - a{Ni{to+)) ^| SVl{C^) 
- SV{Ni) SVl(C^) - A* SVl(CT) SVr{C^) + 0{e) TV(iVi) 
< - SVi(CT)| ( - A* SYr{C^) + 'I^a{Ni{to+))) 
~ SV{Ni) SVl(CT) + 0(e) TV(iVi). 

For any fixed value of t]n G (0, 1), the first term is negative because | SVfl(C^)| < e 
can be assumed a priori smalL The second term is negative and the last one can be 
incorporated into the first term of (|5.6p . This completes the proof of Theorem 1 5. II □ 

6. Quadratic interaction potential (part 2). We continue our investigation 
of interaction functionals adapted to nonclassical shocks. The proposal now is to 
also involve a weight-function in front of the quadratic terms which is chosen to be 
proportional to the difference in (normalized) wave speeds. Such a weighting was 
used successfully by Iguchi and LeFloch [3j as well as by Liu and Yang [11 in their 
analysis of general hyperbolic systems of conservation laws; see also Liu [lOJ for earlier 
pioneering work on systems. 

The interaction functional proposed now for nonclassical entropy solutions shares 
the same monotonicity properties as the functional Qwcak, in the sense that it fails for 
precisely the same interactions. In preparation for the study of systems of conservation 
laws, a complete study of the scalar case requires an examination of the following 
functional. 

So, we introduce 

Qpos{u(t)) ■■= ^6{x,y)<T{x)a{y), 

x<y 

where a{x),a{y) are the generalized wave strengths defined earlier. 



a{x) 



cr(S'), is a wave located at x, 
0, no wave located at x. 



The potential now contains a weight function 9{x,y) defined by 

C* {a{x) ~a{y))^, same monotonicity. 



eix,y) 



otherwise, 



where the constant C, > will be chosen to be sufficiently small and (fc)+ := 
max{fc,0}. More precisely, here we say that two waves have the same monotonicity 
if after mapping all states to the positive region the waves are either both increasing 
or both decreasing. The definition of a is as follows: 
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where 

w+, u+ > 0, 

V>o{ut), u+ < 0. 

This definition allows us to enforce that comparisons only occur between state vari- 
ables belonging to the same region of convexity, i.e. between positive values. Note 
that, in the present formalism, rarefaction are possibly interacting with each other. 

Theorem 6.1 ("Positive interaction" potential for nonclassical shocks). Let 
(f^ be a kinetic function satisfying the properties (Al)-(A3). Consider the weighted 
interaction functional Qpos evaluated for a front-tracking approximate solution u to 
()l.ip . whose discontinuities satisfy (|1.5|) . and with initial data uq G BV(R). // C* is 
taken sufficiently small, then Qpos is strictly decreasing during all interactions except 
in the cases RC-3, CR-4, CC-3 and CN-3 which involve a small wave interacting with 
a crossing shock. 

Proof. We distinguish between several cases of interactions. For the functional 
Q — Qpos we introduce a new form of the decomposition (|4.2p . 

[Q(u(i))] =:[g]i+^a(2/)B(y), 

V 

where the sum over y excludes waves involved in the interaction at time t. Given 
that the waves not involved in the interaction are potentially arbitrary, [Q] will be 
negative if and only if [Q]i and B{y) are negative. As usual, we will rely heavily on 
the inequalities described in the proofs of Theorems 13.11 and 14.21 as well as of Lemma 
1431 

Case RC-1: {R^_^C^)-{C^'). In this case, a single wave is outgoing and we have 
[Q]i < 0. In the proof of Theorem 13. 11 we observed that the states of interest satisfy 

either Ur < ui < it„j, 

or (fioiUr) < Ul < Um- 

For both sets of inequalities, cr(C^) < a{C^') holds and the normalized speeds are 
ordered 

a(C^') < a(C^) <a{Rl). 

Suppose there is a shock at some point y with speed A, then 

B{y) = a {a{C^') - A)+a(C^') - {a{C^) - A)+a(C^) - 

which by the monotonicity just observed must be negative. On the other hand, if a 
rarefaction is located at y, then the estimate on the wave strengths shows that 

B{y) = a{C^') - a{C^) ~ (a(i?|) - A)+(7(i?^) < 0. 

Case RC-2: {Rlci)-{N^' R^'). Our study of this interaction in Theorem B?^ 
provides 

<j{Ni') < <j{Ci), a{R{_') < a{Rl). 



Diminishing functionals for nonclassical entropy solutions 



29 



In turn, it follows that 

[Q]i - a{Ni') <j(R^') - a{Ci) a{Rl) < 0. 

Next, to control B{y) we first consider the relevant speeds. When ipQ{ur) < ui, then 
inequalities (|3.6p are valid and 

a(i?i') < a{Ni) < S(Ci) < a{Rl), 
while when (pQ^Ur) > ui we have 

a(iV^') < a(i?i') < Z{R\) < a{Ci). 
Now, suppose that a shock wave is located at some point y, then 

B{y) =a (S(iVi') - A)+a(7Vi') + a{R[_') 
-a(a(C7i)-A)+a(Ci)-a(4), 

which is negative in view of the above inequalities on the speeds and strengths. A 
similar argument applies if a rarefaction is located at y. 

Case RC-3 : {R\C^)-{nI CT'). This case is more delicate, due to the fact that a 
nonclassical shock is generated from classical waves. The shock strengths satisfy the 
relation (|4.4p which implies 

cj{N^^)+a{C^') <a{C^). (6.1) 

We distinguish between two cases p.7p and (|3.8p depending on the sign of Ur- An 
inspection of the normalized wave speeds associated to the states in p.7[) and (|3.8[) 
demonstrates that 

a{C^') < a{Ni') < a{Rl) and a{C^') < S(C-^) < a(i?^). 

We analyzed an interaction term in Lemma l4.3l similar to 

[Q]i = a(a(iVi') - a(iVi') a{C^') - a(4)r7(C^), 

and we saw that it was possible to take cr(i?^) vanishingly small, while maintaining 

a non-vanishing product cr{Nj. )a{C^'). Therefore, in general we cannot claim that 
[Q]i is negative. 

We now show that in the limit as cr(i?+) 0, B{y) is again positive for some speeds 
A. Assuming a shock located at y, the general form of B is 

B{y) ^C, {a{Ni') - A)+ a(A^i') + (a(CT') - A)^ <j{C^') 
- a {a{Ci) - A)+ <T{Ci) - a{Rl). 

If (f'^{ui) < and inequalities ^3.7\i hold, then take Ur ~^ f^iui) and u„i ui. In 

this limit, a{C'^') < a{C^) < a{Nj_ ) and the quantities a{C^') and (j{Nj_ ) are 

functions only of ui. Therefore taking A e \a{C^),a{N^ )], we find that for those 
limiting values 

B{y) = a {a{Ni') - A)+ a{Ni') > 0. 
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Now, if a rarefaction is located at y, then (|6.ip suffices to show 

B{y) = a(iVi') + a(CT') - a{C^) - (a(i?^) - A)+ < 0. 

Case RN : {RlNi)-{Ni i?i ). As we saw in the proof of Theorem 13. 1[ there are 
two subcases defined by either 

either Lp^ o ip^{ui) < fo o ip^ [um) < ui < u,„, 

or (fil O ip^{ui) <Ui < (f^O ip'{u,n) < Um- 

The normalized speeds satisfy respectively In Case A the speeds satisfy 

a(i?i') < 2(A^i') < a(iVi) < a{R\), or 
a(iVi') < a(i?i') < a{RX) < a(A^4). 

In both subcases, we have cr(-/V± ) < cr{Nj_) and a{R{_ ) < a{Rl_), and the 

inequalities on the speeds a.{Nj_ ) < a{Nj.) and a(i?i ) < a(i?^). 
Using the observations above, we verify that 

[Q]i = a{Ni')a{R^') - a(i?^)a(ivi) < 0. 
If a shock is located at y, then the term 

B{y) =a (a(7Vi') - A)+ a{Ni') + a(i?i') + {a{Ni) - A)+ a(A^i) - a{Rl) 
is negative. When a rarefaction is located at y, then 

B{y) ^a{Ni') + {a{Ri') - A)+ a{R^') - a{Ni) - (a(i?i) - A)+ 
is again negative. 

Case CR-1 : {C^R^)-{Cj. ). Only one wave is outgoing, so we trivially have 
[Q]i < 0. The states of interest satisfy 

ifilium) < Vaiur) < Ul, 
and therefore the normalized speeds are ordered 

a{Ri) <a{Ci) < a(Ci'). 

Recall also that (T(CJj_ ) < cr(CJj_). Assume now that a shock is located at y, then 

Biy) (a(Ci') - A)+ a(Ci') - (a(Ci) - A)+ - 

The least favorable case corresponds to the situation where A = a{C'j.), since only 
one negative term remains when A > a(C£ ) and B{y) is increasing with respect to 
A when A < a(Cj_). The important term in B{y) when A — a(Cj_) is the difference 
in speeds 



f{ui) - f{Ur) f{ui) - f{Urn) 



Ul ~ Ur Ul — 

<C\ur~uJ\ = Ca{R[_). 
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The constant C — C{ui) is bounded above since the total variation of the solution u 
remains bounded. If we therefore assume the a priori bound 

aCTV(w) < 1, (6.2) 

then 

B{y) <C,Cg{r[_)<j{C^^) - <t{rI) ^ -cr(i?i)(l - C,C (T(Ci')) < 0. 

To complete the proof, we assume that a rarefaction is now located at y, and 
immediately see that 

B{y) = a(Ci') - a{Ci) - {a{R[_) - A)+ a{R[_) < 0. 

Case CR-2 : {C\^R\)-{C^ ). The waves are entirely classical and the calculations 
are identical to previous Case CR-1. In particular, we again have to impose 
condition (16.211 on C*. 



Case CR-3 : (Cii?i)-(7Vi'i?i'). In the proof of TheoremO we identified two 
subcases distinguished by the inequalities (|3.10p and (|3.1ip : 

either ip\){um) < V'o ° {'^i) < n < V'o (""»•), 

or (fioiUm) < (PqO (p^{ui) < V5o(Ur) < . 

In both subcases we have 

a{Ni') < a{Ci), (j{r[_') <(j{Ri). 
In the first subcase, the normalized speeds satisfy 

a(Ci) < mm{a{R^_),a{Ni')) < max a(ivi')) <a{R^'), 

while in the second, 

a(i?i) < min(a(Ci),a(i?i')) < max (a(Ci),a(i?i')) <a{N^'). 

The key point is that in both subcases, a(Cj_) < a{N^ ) and a{R^_) < a(i?i ). 
Using the inequalities satisfied by only the strengths, we obtain 

[Q]i = a{Ni')a{Ri') - G(Ci)a{R\_) < 0. 

When a shock is located at y, we can verify that 

B{y) = a (a(iVi') - A)+ a(7Vi') + - (a(Ci) - A)+ a{Ci) - a{Ri). 

The worst case occurs when A = a(C^) < a{N^ ), and then we find that there exists 
a constant C such that 



a{Ni)-a{Ci) 



fiui)-f{ip'oOip\ui)) f{ui)-f{Um) 



Ul — (PqO ip^{ui) Ul — Urn 

< C{^'o ° - "m) = C(a(i?i) - <7{Ri')). 
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We again impose a condition of the form (|6.2p on and find 

B{y) =a {a{Ni')~a{Ci)fa{Ni') - {<j{Ri) - a{Ri')) 
= - (f - Ca cr{Ni')) {<T{Ri) - <T{Ri')) < 0. 
Suppose next that there is a rarefaction at y, then 
B{y) ^a{Ni') + {a{R[_') - A)+ - a{Ci) - (a(i?i) - A)+ 

The largest value of B{y) occurs when A — a(i?i) and then the coefficient of the 
(T(_Ri ) term satisfies 



a{R^') -a{R^) 



Ur - (PqO(P^ {ui ) Ur -Um 

< C{^1 O >^\ui) - U,n) = C {<j{Ci) - <j{Ni')). 



The same arguments as before apply and again allow us to conclude that B{y) < 0. 
Case CR-4 : (Cii?i)-(A^i'cl'). Here, we have 

Now, the normalized speed satisfy 

a{Ri) < a(Cl') < a{Ni'), a(i?i) < a{Ci) < a{Ni') 
and for the strengths we have 

a{Ni')+aiCl') <(T(Ci). 
For the interaction, we obtain a term 

[Q]i^C, {aiNi') -aid'))^ aiNi')aiCl_') - aiC^)aiRi) 



similar to the one analyzed in Lemma 14.31 There we saw that it was possible to let 
(T(i?i) while keeping both (j{N^ ) and cr(Cl) bounded and non- vanishing. 
Therefore, this term is positive for a large class of interactions. 
When a shock is located at y, the general term to look at is 

B{y) ^C,{a{Ni') - A)^a{Ni') +a{a{Cl') - A^aid') 
- C,{a{Ci) - A)+(7(Ci) - f7(i?i). 

Taking again the hmit as cr(-Ri) 0, we find a(Cl ) < a(CJ^) < a{N^ ). Moreover, 

when A = a{Cj.) we obtain expressions (j{Nj_ ) and a{Nj_ ) — a(Cj_) that are 
non-zero and functions only of ui. Therefore 

B{y) = C,{a{Ni')^A)+a{Ni') > 0. 

When a rarefaction is located at y, then 

B{y) = a{Ni') + a{Cl') - a{Ci) - C,{a{R^) - A)^a{R^) < 0. 
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Case CC-1 : (Cj_C^)-(C^'). There are two subcases determined by the sign of Ur : 

(foiUr) <Um < Ui, Or Ur < Um < Ui- 

In both subcases, we have the identities 

cj{C^') = a(Ci) + a{C^), 
a(C^) < a{C^') < a{C\.). 

Only one wave is outgoing so [Q]i < 0. 

Suppose a shock is located at y, then the largest value of 

B{y) = a{a{C^') - A)+a(C^') - C\{a(Ci) - A)+ct(C|) - C»(a(C^) - A)^a{C^) 
occurs in the interval A e [o(C^), o(C^')], leaving us with 

B{y) = a{a{C^') - A)+<7(Ci') - a{a{Ci) - A)+c7(C|). 
In the subcase where Ur > 0, the definition of a and a allows us to find 

B{y) =a (^{f{ui) - fiUr)) - A{UI - Ur) - (fiui) - f{Um)) + A(«i - Um)) 

=C*(a(C^) - A)ct(C^) < 0. 

If Ur < 0, then similar computations show that B{y) is still negative. 
When a rarefaction is located at y, then 

B{y) = a{C^')-a{Ci)-a{C^) = 0. 

Case CC-2 : (CicT)-(C^'). Two subcases appear depending on the sign of Ur- We 
will assume throughout that Ur > 0, leaving us with the states 

Ur < (fioiUm) < Ui, 

and therefore the two identities 

a{C^) = a{Ci) + a{C% 
a{C^) <a{C^') <a{Ci). 

We recognize immediately the same identities as in the Case CC-1. As quick 
verification shows that the same arguments demonstrate that [Q]i < and [Q]2 < 0. 

Case CC-3 : {C\_C^)-{Ni C^' ). We can distinguish between two subcases : 
either y>o("r) < < '/'o ° V^i^i) < ui, 

or foiUr) < ° <Um< Ul- 

In both subcases, one obtains 

a(C^) < a{Ni') and a{C^') < a{C\_). (6.3) 
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Unfortunately, the conditions above on the shock speeds do not preclude the 
possibility that a(C^) < a{Nj_ ) which would imply that 

B{y) ^C,{a{Ni') - Afa{Ni') + a(a(CT') - A)+a(Cf ') 
- C,{a{Ci) - A)+a(C|) - C,(a(C^) - A)^a{C^) 

was positive. On the other hand, because [V] = we have that B{y) = when a 
rarefaction is located at y. 

The sign of [Q]i is also less than helpful. Using the ideas seen in Lemma [4.31 for the 
same term, one notices that when Um equals the transitional state tp^ o (p^{ui) 
between the two subcases, then the change 

[Q]i = C,{a{Ni') -a{C^'))+aiNi')aiC^') - (a(Ci) - a(C^))+(7(Ci)a(C^) 

vanishes because the quadratic terms and the coefficients are equal. It is clear that 
for some set of waves close to this transitional value, [Q]i will also be positive. 

Case CN-1 : {C\_Nj_)-{Cj. ). The states are ip^ o Lp^{um) < Um < ui and these 
provide us with the basic inequalities 

aiCi')^cT{Ci)+aiNi), 

a{Ni) <a{Ci') <a{Ci). 

Under these conditions, the analysis is similar to the CC-1 case. 

Case CN-2 : (CiiV^)-(C|'). The important states are ordered 
(p^{u,n) < foi'^m) < ui and therefore 

a{Ci') = a{Ci)+a{Nl), 
a(Nl) <a{Ci') <a{Ci). 
Again, the analysis is similar to the CC-1 case. 

Case CN-3 : {cIn\_)-{nI' Cl'). There are two subcases defined by the ordered 

set of states 

either i^q ° ^^(um) < < V'o ° < '^h 

or O ip^{Urn) < ¥^0 ° M < < '>^h 

where the second case corresponds to a weak incoming shock C^. It is easy to verify 
that the shock strengths verify 

a{Nl) < a{Nl') and cr{C^') < cr(Ci). (6.4) 

In both subcases, one obtains the unfortunate inequalities 

a{Ni) < a{Ni') and a{C^') < a{Ci). (6.5) 

The analogue to 

[Q]i = C4a{Ni')-a{Cl'))^<j{Ni')a{Cl_') - C^aiCi) - a(A^i))+a(Ci)a(A^i) 
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for Qwcak was analyzed in Lemma 14.31 and seen to be positive for sufficiently strong 
waves C\_. For the functional Qpos, [Q]i = when u„j = ° V^i'^i) because at that 
critical value, all inequalities in (|6.4p and (|6.5p are equalities. Without going into 
the details, it easy to verify that this implies that for a non-empty set of waves, [Q]i 
will be positive. 

Suppose a shock is located at y, then 

B{y) ^a{a{Ni') ~ A)+a(7Vi') + (a(Cl') - A)+a(Cl') 
- a(a(Ci) - A)+a(Cii) - C^aiNi) A)+c7(7Vi). 

When we are in the first subcase (a strong C\_), the change B{y) is clearly positive if 
AG[a(Ci),a(iVi')]. 

Finally, suppose a rarefaction is located at y, then 

Biy) = a{Ni') + ct(CI') - a(Ci) - a{Ni) = 0. 

Case NC : {N^C^)-{C^'). The defining states are Ur < (fQ o ip^{ui) < ui and we 
have 

a{C^) <a(C^) <a{Ni). 

See the CC-1 case. 

Case NN : (7V_^iV^)-(C|'). The analysis of the NC case also apphes here since NN 
interactions are limiting situations obtained as Ur — s- 
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